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NOTATION 


a 

h 

C 

Cx 

c 

D 

D 1  +  iD  2 
e(0) 

F 

f(k,a,z) 

G(k,) 

9 

g{k) 

g'{k) 

H(m) 

»2 

«3 

«4 

h 

l{r),w) 


sinh  Ka  cos  K a 
cosh  Ka  sin  Ka 
cosh  Ka  cos  Ka 
sinh  Ka  sin  Ka 

Height  of  hydrofoil  above  sea  bed 

Semispan  of  hydrofoil 

Average  lift  coefficient  of  hydrofoil  * 

(^2D,  •-/V,£>2)/  (0,2  +  0/)  • 

Chord  of  hydrofoil 

Subscript  to  designate  line  doublet 

Complex  denominator  of  Wj  («>,»;)  * 

Functions  of  0  occurring  in  K  (ij)  • 

Froude  number  based  on  water  depth 

Functions  of  a,  and  z  occurring  in  K  (17) 

Undetermined  pu^ameter 

Acceleration  of  gravity  » 

sinh  kh  cosh  kh 

kh  cosh  2  kh 

Undetermined  parameter 
sinh  2Kh  cos  2Kh 
cosh  2Kh  sin  2Kh 
cosh  2Kh  cos  2Kh 
sinh  2Kh  sin  2Kh 
Average  water  depth 

Imaginary  part  of  kernel  function  K{j\)  * 
Real  part  of  kernel  function  K  (q)  * 


•  Primed  qnantities  ere  used  for  the  eopercriticel  rsnge  of  Froude  number  and  unprlmed  quantities  for  the 
•ubcritical  ranges 


K  Integration  variable 

^  W)  Kernel  of  an  integral  defining  a  potential  function  • 

k  Integration  variable 

^0  Pole  on  the  real  /(--axis 

LijO  Kernel  of  an  integral  defining  a  potential  functio'h  * 

Coipplex  kernel  function 
m  Integration  variable 

A/j  +  Complex  terms  in  the  numerator  of  itfj  (w,e) 

p  Pressure  change  in  the  fluid 

Q  Unity  if  F  <  1  and  (-  1)"  if  F  >  1 

+  Sj 

/?2  x'^  +  ^2 

T 

Sj  (ti  -  +  (2 nh  +  z  -  a)^ 

^2  (i;  -  y)^  +  (2nA  +  z  +  a)^ 

T  (i;)  Kernel  of  an’integral  defining  a  potential  function 

U  Free-stream  velocity 

V  Subscript  to  designate  semi-infinite  vortex  sheet 

Either  or  w_ 

®  cos  0  +  (rf-y)  sin  6 
X  cos  0  -  (rj-y)  sin  0 
X  Longitudinal  coordinate 

y  Horizontal  lateral  coordinate 

sinh  K z  cos  l(z 
^2  cosh  Kz  sin  Kz 

^  t  • 

^3  j  cosh  K z  cos  Kz 

^4  sinh  Kz  sin  Kz 

z  Vertical  coordinate 

r'('?)  Spanwise  distribution  of  vorticity  , 

C  Wave  elevation  ^ 

7  distance  measured  along  the  span 


IV 


Shape  parameter 

Horizontal  doublet  distribution  21ong  the 
Distance  in  the  ^'direction 
Density  of  fluid 

Total  potential  function  ^ 

Perturbation  ^tential  function  * 


ABSTRACT 


This  report  shows  that  singularities  for  a  semi-infinite  vortex  sheet  , 

a  horizontal  line  doublet,  both  of  finite  span,  may  be  used  to  represent  the  vor- 
ticity  and  displacement  produced  by  a  finite  span  hydrofoil.  The  boundary  con¬ 
ditions  which  must  be  satisfied  when  these  singularities  move  in  shallow  water 
are  derived.  Expressions  for  tife  potentiaj  functions  are  obtained  for  both  the 
subcritical  and  supercritical  ranges  of  Froude ^number.  Explicit  expressions  for® 
the  X-,  y-,  and  ^-derivatives  of  the  potential  functions  are  also  presented.  Meth¬ 
ods  for  determining  the  strengths  of  the  singularities  are  discussed. 

INTRODUCTION 

As  an  extension  of  Wu’s  theoretical  study  of  finite  span  hydrofoils  mgving  beneath  the 
surface  in  water  of  infinite  depth, ^  expressions  are  derived  here  for  the  potential  functions 
associated  with  finite  span  hydrofoils  in  water  of  finite  depth.  These  expressions  may  be 
obtained  for  any  span,  water  depth,  or  depth  of  submergence,  and  for  any  velocity  in  the  sub- 
critical  and  supercritical  ranges  of  Froude  number.  In  additiob  to  the  potential  functions, 
explicit  expressions  are  derived  for  the  x-,  y-,  and  a-derivatives  of  these  functions.  With 
the  results  of  Wu’s  report  and  those  of  otlrer  investigators,^’^  these  functions  may  be  used 
to  determine  wave  profiles,  wave  resistance,  velocity  components,  pressure,  and  other  pa¬ 
rameters.  , 

A  few  theoretical  studies  h&\e  been  made  of  the  wavemaking  produced  by  singular¬ 
ities  moving  beneath  a  free  surface  in  shallow  water.  Lunde^  outlined  the  method  for  find¬ 
ing  the  wavemaking  of  a  source  in  shallow  water,  .nd  Pond®  and  DiDonato®  have  made  a 
complete  analysis  of  the  problem  for  t  .e  subcritical  range  of  Froude  number.  A  further  dis¬ 
cussion  of  problems  of  this  typie,  including  a  finite  span  hydrofoil,  is  given  in  Reference  7. 
None  of  these  works  gives  a  complete  analysis  for  the  supercritical  range  of  Froude  number. 

When  a  singularity  moves  horizontally  between  a  rigid  bottom  and  a  free  surface,  the 
potential  function  is  very  sensitive  to  the  value  of  the  Froude  number  based  on  water  depth. 
As  the  velocity  increases,  the  wavemaking  builds  up  and  becomes  very  severe  when  the 
Froude  number  approaches  unity.  When  the  Froude  number  exceeds  unity  the  singularity  is 
moving  faster  than  the  wave  velocity  and  wavemaking  effects  diminish.  This  phenomenon  is 
common  where  ships  enter  shallow  water  and  the  resistance  suddenly  decreases  when  the 
velocity  exceeds  the  critical  speed.  Although  most  surface  ships  seldom  exceed  critical 
speed  in  normal  operation,  this  range  may  be  important  in  the  operation  of  hydrofoil  boats. 

Wu  has  shown  that  the  singularity  which  produces  the  same  vorticity  as  a  hydrofoil 
is  a  s'emi-infinite  vortex  sheet  of  finite  span,  extending  from  the  position  of  the  hydrofoil 

^References  are  listed  on  page  29.  ^ 
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to  an  infinite  distance  downstream.  In  addiUon  to  the  vorticity,  the  hydrofoil  produces  a  dis-® 
placement  effect.  At  some  distance  from  the  hydrofoil,  a  •horizontal  line  doublet  of  finite  span, 
whose  strength  is  proportional  to  the  water  displaced,  is  adequate  to  represent  the  displace# 
ment  effect.  In  shallow  water  both  singularities  must  satisfy  boundary  conditions  at  the  sea 
bed  and  at  the  free  surface.  Furthermore.^e  water  surface  far  ahead  of  the  hydrofoil  is  not 
disturbed.  In  this  report  expressions  are  derived  for  the  potential  functions  of  each  of  the 
singularities.  Because  the  boundary  conditions  and  methods  of  analysis  are  the  same  for 
each  casd,  general  considerations  will  be  discussed  first.  Finally,  a  method  for  determining  ® 
the  strengths  of  the  Angularities  will  be  discussed. 

BOUNDARY  CONDITIONS  FOR  THE  THEORETICAL  ANALYSIS 

• 

The  hydrofoil  may  be  considered  fixed  in  an  inertial  system  in  which  there  is  a  uniform 
velocity  f/  at  z  =  -  ~  flowing  in  the  positive  avdirection.  As  shown  in  Figure  1,  the  hydro¬ 
foil  with  chord  c  and  semispan  6  is  fixed  at  *  -  0,  s  =  a  with  its  span  extending  from  y  =  -  6 

to  y  -  6.  The  origin  of  the  coordinate  system  is  at  the  sea  bed  and  the  plane  a  =  A  is  the  un¬ 
disturbed  water  surface.  *  * 


Figure  1  -  Coordinate  System  for  a  Hydrofoil  in  Shallow  Water'  ^ 

The  potential  flow  about  a  hydrofoil  may  be  obtained  by  representing  the  hydrofoil  by 
several  singularities  extending  along  the  span.  The  steady-state  potential  function  is 


2 


4>  (*,  y,  b)  ~  U  X  +  <f>  (ar,*y,  2) 


o 


where  <f>  is  the  perturbation  potential.  If  disturbances  in  the  flow  are  small,  interaction  terms 
may  be  neglected,  and  is  the  sum  of  the  potential  functions  for  each  of  the  singularities 
used^o  represent  the  hydrofoil.  The  change  in  pressure  p  produced  at  any  point  in  the  fluid 
fs  obtained*from  Bernoulli’s  equation.  Along  any  streamline  *  ^ 

1 

p  {x,  y,  b)  =  -  p  U  - - pg  {s-h)  -  —  p  (grad  [2] 

dx  ^2 

» 

where  p  is  the  density  of  the  fluid, 

g  is  the  acceleration  of  gravity,  and 

B  is  the  position  of  the  streamline  in  the  undisturbed  flow. 

In  the  linearized  theory,  the  last  term  in  the  Eqi^tion  [2]  may  be  neglected.  Along  the  stream¬ 
line  represented  by  the  free  surface,  the  pressure  change  produced  is  zero  and  the  wave  el¬ 
evation  is  given  by 


9  z=h  • 


If  the  sea  bed  is  considered  as  a  rigid  surface,  there  is  no  flow  normal  to  this  surface 
and  the  first  boundary  condition  is 


Boundary  Condition  I 


© 


For  disturbances  at  the  free  surface  which  are  small  compared  with  the  wave  length,  the 
normal  velocity  of  the  surface  is  the  same  as  the  vertical  velocity  of  the  fluid  particles. 
Then  • 


.0 

\  dx  dz  /  , 


If  Equations  [3]  and  [5]  are  combined,  the  second  boundary  condition  which  applies  at  the 
free  surface  is  obtained; 


Boundary  Condition  II  (  — ^  +  F^h  — 

\  ds  dx^  / 


dx  /  z  =  h 


where  F  is  the  Froude  number  based  on  water  depth 

•  U 

F  = 

•  V  fffi 


% 


3 


o 


As  the  free  surface  far  ahead  of  the  hydrofoil  is  undisturbed,  the  wave  elevation  and  hence 
the  potential  itself  must  vanish  far  upstream.  Therefore  the  third  boundary  condition  is* 


Lioundary  Condition  III 


for  a:-*-oo,  ^  =  0 


In  addition,  Ih,  polenU.I  and  it,  deriv.tivo.  are  finite  ev»ywh„e  in  the  region  except  at  the 
position  of  the  singularity.  ^ 

If  the  hydrofoil  is  represented  as  a  lifting  line  with  finite  displacement  and  its  vortex 
wake,  two  singularities  are  required  to  represent  the  foil  in  an  unbounded  Huid.  The  circu¬ 
lation  and  trailing  vortex  sheet  are  represented  by  a  semi-infinite  vortex  sheet  of  finite  span. 

he  displacement  of  the  foil  ,s  represented  by  a  horizontal  line  doublet.  As  the  problem  is 
linear,  the  total  potential  is  the  sum  of  the  potentials  for  the  two  singularities.  Both  poten- 
t.als  must  satisfy  the  three  boundary  conditions  independently.  The  solution  for  each  poten¬ 
tial  IS  built  up  from  the  potential  of  that  singularity  in  an  unbounded  fluid.  The  function 
IS  the  potential  of  the  singularity  in  an  unbounded  fluid  plus  iti  image  in  the  sea  bed  such  ' 
that  Boundary  Conditions  I  and  III  are  satisfied.  The  function  is  determined  in  such  a 
way  that  the  potential  satisfies  Boundary  Condition  II  as  well  as  I.  The  potential 

<t>,  +  IS  found  to  consist  essentially  of  the  sum  of  three  terms  which  will  be  designated 
as  ^  and  ,^3.  In  the  subcritical  range  of  Froude  number,  ,  is  the  total  potential  for  * 

zero  Froude  number,  4,^  is  a  short-range  potential  which  decays  rapidly  to  zero  a  short  dts- 
ance  rom  the  singularity,  and  4,^  is  the  potential  associated  with  wavemaking.  In  the  super¬ 
critical  range  of  Froude  number,  4>,  is  the  total  potential  for  infinite  Froude  number  but  the 
other  terms  have  the  same  significance.  Boundary  Condition  III  i‘s  needed  to  exclude  wave¬ 
making  far  upstream  from  the  hydrofoil  and  to  determine  the  explicit  form  of  4>^  for  positive 

•and  negative  values  of  ..  When  all  boundary  conditions  are  satisfied,  the  total  potential 
function  IS  «  ^ 

^  4>i  + =  4>i  +  4>2  + 4>3  •  [9] 

POTENTIAL  OF  A  SEMI-INFINITE  VORTEX  SHEET  OF  FINITE  SPAN 

The  potential  of  a  semi-infinite  vortex  sheet  at  a  =  a  and  its  image  at  g  =  -  a  is 
obtained  from  Wu’s  report*  * 


f  |-^  _  a,  -| 

ib  (7-y)^ +(2-0)2  L  +  (,-y)2  +(g_a)2]V4j 


2  +  a  r  »  r  (ij)  dif 

-b  (’?-y)2  +  (g-a)2 


[x^  +  (»?-y)2 


+  (2  +  a)^]^  j 
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van- 


where  r(T;)  is  the  distribution  of  circulation  along  the  span.  It  is  easy  to  see  that  <f>^ 
ishes  for  a;  =  -  ~  and  its  a-derivative  vanis.hes  at  a  =  0.  This  ^nation  also  has  an  integral 

form  which  is  derived  in  Appendix  A  and  which  is  the  real  part  of  • 

*  • 

I  ^  #****  • 

4^  Jc*'"<’?-r)[sgn  (a-o)  rfm 

^  L  ''  -A  0 


-b  “„/2  "  “  ' 


where 


to  =  X  cos  0  +  {ri-y)  sin  0 


*  ^ 

and  sgn  is  the  signum  function  whose  value  is  +  1  when  the  argument  is  positive,  —1  when 

the  argument  is  negative,  and  0  when  the  argument  is  zero.  The  line  Uirough  the  integral  sign 
indicates  thr^t  the  principal  value  is  to  be  taken.  In  this  case  the  integral  must  be  evaluated 
from  -  0/2  +  <  to  n’/2  -  t ,  where  e  approaches  zero  in  like  manner  at  the  two  limits. 

To  satis^  Boundary  Condition  II  at  the  free  surface,  an  undetermined  potential  function 
added  to  This  function  wfll  be  given  th^e  same  x  and  y  dependence  as 

and  must  satisfy  Boundary  Condition  I.  If  <^y  is  written  in  terms  of  undetermined  parameters 
H(m)  and  G(k,6),  then  the  form  for  a  >  o  is  * 


He  —  J  l'{ri)dri  J e^"' sinh  mo  te"'"''  +  //(m)  cosh  ttjz]  (foi 

+  -7  /  - -  sink  kale-'^  +  G(k,0)  cosh  ks]  d/c>  [13] 

-n/2  0  J 

.  A 

C 

The  integrands  of^each  of  these  integrals  must  vanish  when  they  are  substituted  into 
Equation  [6],  the  expression  for  Boundary  Condition  II.  Then 


//(m)  = 


sinh  mh 


G  {k,  (9)  =  - 


II  +  M  cos^  e] _ 

cosh  kh  [kh  cos^  6  -  tanh  kh] 
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When  the  Froude  number  is  zero,  G(k,  0)  becomes 


Goik,  d)=—~ 
sinh  kh 


The  potential  for  zero  Froude  number  is 


1  r  1  ? 


(r/)]  rfr/  [17] 


where  A  is  introduced  to  make  the  kernel  functions  dimensionless  in  length.  The  functions 
('l)  ^nd  L  y  ^  (r/)  have  the  dimensionless  forms 

oo 

Ty(T,)^  Re—  y'c*'"  ('7-y>rsgn  (2-a)  -“I  -  c-'"  (» +“)  ,  2e~'"*  sinh  mg  cosh  wal, 

0  ^  sinh  mh  J 


V.  y 


h  do 


J  cos  0  J 
—!>/ 2  0 


[o  “A  H 

sgn  (z-a)  ""I  _  g-*  (i  +  o)  ^  g  sinh  ka  cosh 

sinh  kh  J 

If  the  relations  given  in  Appendix  A  are  used,  these  expressions  have  the  nonintogral  forms 


.<n).i  y  f  — 

“  L ('/-y)^  +  (3+2 


3  +  2  nA  +  a 


nh-af-  {rj-yf  +  (B  +  2nh  +  af 


^  V.  (7)  =  — 


\  Q  ^  ^ _ {z  Jr2nh-  a)  x _ 

n  =  -oxi  [(7~y)  +  (3  +  27iA-a)^]  [x^  +  (rj-y)^  +  (z  +  2nA-a)^]^ 

(3  +  2nA  +  a)  ar  ^ 

[(7-y)^  +  (g  +  27iA  +  a)2]  [ar2  +  ^  (2  +  2nA  +  a) 2]'^  [■ 


where  Q  is  unity.  For  intermediate  Froude  numbers  in  the  subcritical 
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[22] 


G  (A,  6)  =  Go  (A,  0)  - 


cos'^  0 

sinh  kh  cosh  kh  {khF^  cos^  d  -  tanh  kh\ 


•  1  r 

^Vj  ^Vjj  "  ^  7  ^ 


where 


^  n/  A  «» 

^  A  /•  /• 

Ky  (ri)  =  -  Re  - 7  7”  ~7 


e**""  sinh  A:a  cosh  Arg  kdk 


cos  0(f0  -j- - [24] 

7  sinh  A:A  cosh  kh  [khF^  cos^  d  —  tanh  M] 


For  the  hypothetical  case  of  infinite  Froude  number 


{k,  0) - 

e  “  '  ’  '  cosh  kh 


b  b  ^  . 

-6  *  -fc 

where  Ty  (ri)  is  already  defined,  and  Ly^  (rj)  is  the  same  as  Ly^  (rj)  exctpt  for  an  alternating 
sign  in  the  infinite  sum.  Therefore,  Ly^(rf)  is  given  by  Equation  [21]  if  Q  is  given  the  value 

<?  =  (-l)'’  •  *.  [27] 

For  intermediate  Froude  numbers  in  the  supercritical  range 


G  (/c,^0)  =  (k,  0)  + 


1  kh  F'^  cos^  0 

-  1 - 

sinh  M  cosh  AAL  M cos^  <?  *  tanh  M 


and  the  potential  is 


where 


fir  <• 


(’7)  =  (v)  +  A"  (17) 


7 


In  this  expression 


n/7.  *00  .• 

t  ’  i-\  D.  ^  C  C  «  sinh*acosh/t2rf/t 

7i  J  ;^1,J  -  .i„hM,o.hu  [311 

-n/2  0 

Therefore  Ly  (1,)  is  twice  the  sum  of  the  terms  of  Ly  (77)  for  which  (/=-!. 

%  • 

Since  the  potentials,  associated  with  the  free  surface  condition,  for  the  two  ranges  of 
Froude  number  are  similar  in  form,  the  following  discussion  will  apply  to  both  functions.  The 
range  of  integration  over  negative  values  of  6  may  be  changed  to  positive  values  if  w  is  re¬ 
defined  as 


=  X  cos  0  +  (77-y)  sin  0 
=  X  cos  0  -  (t)  -  y)  sin  0 


Then  K  y  (rj)  becomes 


[32] 


17/2 


A'j,  (77)  ^  He  i  J  e  {0)  [W  {ic^,  6)  +  M  {w_,  e)]  dd 


[33] 


where  ,W  (w,  0)  is  the  complex  integral 


M  (w,  0) 


09 

=  7/ 


/(/',<2,2)  dk 


0  g  {k)  cos2  6 


tanh  kh 
kh 


[34] 


In  this  equation  to  is  written  for  either  or  w_.  The  functions  f  {k,a,z),  e  (0),  and  g  (/t) 


are 


f  {k,a,z)  =  sinh  ka  cosh  kz 
e  (6)  =  cos  6 
g  {k)  =  sinh  kh  cosh  kh 


[35] 


The  integrand  of  M  {w,  0)  has  a  simple  pole  when  k  =  where  is  defined  by  the 


equation 


k^h  F^  cos^  6  =  tanh  k^h 


[36] 


This  integral  may  be  evaluated  by  means  of  a  contour  integration  in  the  complex  /t-plane. 
The  method  is  described  in  Appendix  B. 
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The  contour  integration  yields  a  value  of  K y  (rf)  which  is  the  sum  of  two  terms,  K, 


and  In  the  subcritical  range  of  Froude  number 


f<  y  (n)  =  /a  +  >2  (V.’O  *  [37] 

•  ^ 

# 

where  /j  (i?,'/')  is  defined  in  Equation  [15B],  \ppendix  B,  in  which 


'5) 


/Vj  -  AyZ^  ^2^4 
/Vj  =  AyZ^  +  A^Z  ^ 
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[38] 


The  values  of  the  terms  used  in  Ny  and  /Vj  are  given  in  Equations  [17B]  and  [18B].  In  the 
supercritical  range  of  Froude  number  the  integrals  defining  h' y^  (q)  have  singularities  which 
occur  when  F  cos  0  =  1.  Therefore  it  is  necessary  to  find  the  Cauchy  principal  value  at  the  • 
singularity. 

K  y^  (r;)  =  L  y^  (rf)  +  R  y^  (rj)  [39] 


It  might  appear  to  be  more  natural  to  combine  Ly^  (q)  with  K  y^  (r,)  to  give  Ky^  (q),  the  zero 
Froude  number  value.  However,  from  the  few  experimental  and  theoretical  studies  which  have 
been  made  with  hydrofoils  operating  at  Froude  numbers  slightly  greater  than  unity,  the  term* 
('?)»  as  given  here,  is  a  good  approximation  for  the  flow  near  the  hydrofoil.®  The  approx¬ 
imation  is  even  better  for  higher  Froude  numbers.  This  indicates  that  the  other  nearfield 
» 

terms  which  are  contained  in  h  y^  (q)  are  of  secondary  importance  and  probably  may  be  ne¬ 
glected  in  many  computations. 

The  other  integral  terms  obtained  in  the  contour  integration  are  of  the  form  (q,^). 
Equations  [22B]  and  [26B]  in  Appendix  B.  As  these  forms  depend  upon  regions  in  which 
and  w_  are  positive,  they  will  be  written  explicitly  for  positive  and  negative  values  of  x, 
using  the  resiflts  gi>ten  in  Table  1  of  Appendix  B.  In  the  subcritical  range  of  Froude  number 
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o 


0 


ar>0,  7;-y>0 


L7"  »r/2-6lj 

I  g  (0)  cos  k^w^dd  +  /  ^  (0)  cos  kQW_dd  I 

.•  n  ®  J 


x<0,  g  —  y<0 


r  (0)  cos 

;r72-0, 


g~y 

where  tan  0,  =  | - ,  and 

X 


<7(0)  = 


k^h.  cos  0  sinh  k^a  cosh  k^z 
sinh  k^h  cosh  k^h  cos^  0  -  sech^  AqA] 


In  the  supercritical  range  of  Froude  number 


x>0,  g-y>0 


V.  (»?)  =  -  ^ 


/  /(^o)  k^w^dkQ  +  /  /  (Aq)  cos  Apt-WAo  [43] 

-J 


X<0,g-y>0 


oo 

=  J  /(^o)  cos 


where  A is  the  value  of  Aq  for  which  0  =  77/2  -  0j  and 


M*o)  = 


AqA  sinh  AqO  cosh  AqZ 
sinh  AqA  cosh  AqA  sin  0 


For  negative  values  of  77  -  y,  and  are  interchanged. 

The  potential  <f>y^  which  contains  the  kernel  Ky  (rj),  is  a  nearfield  term  that  decays 

rapidly  with  the  distance  from  the  singularity.  The  potential  <f>y  ,  which  contains  the 

®  3 
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kernel  Ky^  (r,),  is  a  wavemakinp  term  which  accounts  for  the  waves  downstream  from  the  hydro¬ 
foil.  The  waves  extending  upstream  from  the  hydrofoil  are  of  high  frequency  and  Secay  rapidly 
with  distance.  Finally,  the  total  potential  function  for  the  semi-infinite  vortex  sheet  in  shal¬ 
low  water  is 

® 

^  6 

^  y*  *  (*1)  ('?)  +  (’?)  +  (’i)l  [46] 

®  -6 

The  X-,  y-,  and  ^-derivatives  of  are  readily‘obtained  from  the  derivatives  of  the 
kernel  functions.  The  three  derivatives  of  are  obtainedTfrom  Equations  [20]  and  [21]. 

If  ffj,  R ^11  and  Sj  are  introduced  for 

•  •  • 

+  (V~y)^  +  (2  +  2nh-a)^ 

0  2  2  9  9  t^7] 

^2  =  -  yr  +  (z  +  2nh+  a)^ 

+  (2  +  2nA-o)^ 

o2  2  9 

^2»“  i^~y)  +  (3  +  2nA  +  o)^ 


the  derivatives  ot  T y  +  L  y  are 


d  {Ty  ^Ly  )  ~ 


nh  -  a  z  +  2nh  + 


a 

—  [49]  . 


^ {Ty  +  Ly  ) 

i  •  A 


{i-y) 


Z 


f  2  (3  +  2nA-  a)  [ 

’l  .  *  1 

1  1 

1  +  -  1 

ri.ifi 

Q‘x  (2  +  2n4-a) 

+ _ 1 

('*  (z  +  2  nA  +  o) ^ 

L  R,  J 

S  2  j 

* 

11 

• 

« 

•  • 

% 

d  {Ty  +  Ly  )  2 

•a  =  * 

-  (2  +  2nA-o)2  r  •Qx'] 

dz  2 

1  S'  L  *.J 

* 

n  *  —00  1 

(77-y)^  -  (2  +  2nA+ 

Qx  (z  +  2nh-a)^  Qx(z  +  2nh  +  a)^ 

L  «2J 

S2/?3  siRi 

The  derivatives  of  the  other  kernel  functions  are  obtained  from  the  derivatives  of  Ky  (77)  in 
Equation  [33].  The  ^-derivative  is  the  real  part  of  ' 

3Ky  ‘ 

A  =-•/('?, »"+)  J  ’  e  (0)  [V  0)  +  /W  (M_,  0)1  dd  [52] 

0* 

•  • 

where  e  (0)  is  the  product  of  cos  0  and  its  former  value.  The  y-derivative  is  • 
dKy  * 

^  =  J  -  J  J  e  {e)[M  {w^,d)  -  M  {w_,e)]  dd  [53] 

.  0 

^  • 

whe^e  e(0)  is  the  product  of  sin  0  and  its  former  value.  In  both  cases  the  new  value  of 
f{k,a,z)  is  •  *  • 

,  f(k,a,z)  =  kh  sinh  ka  cosh  kz  [54] 

Expressions  for  the  functions  (17, ic)  and  may  be  found  in  Appendix  B.  The  func¬ 

tions  N  J  and  A'j  needed  for  cfj  (’I* '‘f®  obtained  by  replacing  k  hy  K  (1  +  i)  in  Equation  [54]. 
Then  the  real  and  imaginary  parts  are 

=  A'A  [AjZj  -  A2Z4  -  -  A2Z3] 

=  Kh[A^Z^- +  A^Z^] 


The  2-derivative  of  K y(ri)  is 


A  =  /  (77,  tC^)  +  /  (77,  M7  J 


The  function  e(0)  is  unchanged  but  f{k,a,z)  becomes 

f{k,a,z)  =  kh  sinh  ka  sinh  kz 
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^xpre|sions  for  (r/,  w)  and  w)  may  be  found  in  Appendix  B.  The  functions  and 


required  for  (ri,  vu)  are 


^J  -  O  [i4-,Z,  -  AjZj  -  AjZj  -  ^2^1^  ^ 

'  A'j  =  A'A  [AjZj  -  AjZj  +  AjZj  + 

Explicit  forms  for  the  three  derivatives  of  Ky^{r))  for  positive  and  negative  values  of  x  are 
readily  obtained  from  the  derivatives  of  the  functions  given  in  Equations  [40]  through  [45]. 

POTENTIAL  OF  A  HORIZONTAL  LINE  DOUBLET  OF  FINITE  SPAN 

.  ♦  ' 

The  potential  of  a  horizontal  line  doublet  of  finite  span  at  3  =  a  and  its  image  in  the 
sea  bed  at  2  =  -  a  is 


,  - 


^  =  xA  j  fi(r)) 

-b 


[x^  +  (v-y)^  +  [x^  +  (rj-y)^  +  (2  +  a)2]3/2 


drj  [59] 


where  h(ti)  is  the  doublet  strength  per  unit  length.  This  potential  vanishes  for  large  values 

of  X  and  its  2-derivative  vanishes  at  the  sea  bed.  The  integral  form  for  is 

•  *  / 

« 

^  b  n/2  oo 

=  J  il)  dv  J  cos  Ode  J  e‘*“'  [e“*l*“°l  c“*^*'*'“>]  AdA  [60] 


where  w  is  defined  in. Equation  [12]. 

To  satisfy  the  free  surface  boundary  condition  of  Equation  [6],  it  is  necessary  to  add 
a  potential  which 'contains  an  undetermined  parameter  G  (A,  0).  As  in  the  case  of  the 

senfi-infinite  vortex  sjieet,  G  (A,  0)  must  have  the  same  x  and  y  dependence  as  <f>jy  and  must 
satisfy  Boundary  Condition  I.  For  points  near  the  surface  where  2  is  greater  than  a,  the 
potential  is 

2A  * 

9^0/ +  J  h{ri)dr,  J  cos  OdO  J  c‘*“'  cosh  Aa  [c-**  +  f;((9,  A)  cosh  A2I  ArfA 
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When  Boundary  Condition  II  is  applied,  G{k,  6)  has  the  same  form  it  had  for  the  case  of  the 
semiynfinite  vortex  sheet,  Equations  [15]  through  [^8].  It  is  therefore  possible  to  write  down 
the  solution  at  once.  • 


2 

"  °  A  /  D  D  Sn)  *  D 

•^6  ^ 

O 

In  this  equation  the  values  of  (,)  for  the  two  ranges  of  Froude  number 

•  1 


[62] 


are 


^2  ^ 1 2  oo  .  • 

j  ('?)  =  —  f  cos  OdO  f  c'**"  r  I*  “  “I  +  e“*(*  ^  t  2e  cosh  ka  cosh 

:l„/2  0  ^  ‘  sinh  k/i  J 


[63] 


,2 


^  n  —  f  cos  OdO  f  e  **“'[' +  ^  •  cosh  ka  cosh  Aa  "j 

‘  .  2"*  4/2*  {  - J 


kdk 


[64] 


_  These'functions  have  the  nonintegral  form 


—  ^  Q  j  - - 

^  1  +  (»?-y)^  +  (z  + 


27iA-a)2]3/2  [3.2  ^  ^  (2  +  2nA  +  a)2]3/2 

« 

[65] 


case* 


where  Q  is  unity  in  the  subcritical  case  but  has  the  value  (-I)"  in  the  supercritical 

The  functions  K  jj  ^{rj)  and  A'q^Ctj)  are  obtained  from  the  real  part  of  the  double  integral 


n/  2 

(n)  =  /  (tc^,  r,)  +  /  («,_,  r,)  i-J  e  (0)  IM  0)  +  Id  (w_,  0)]  d6 


do  [66] 


where  M  {w,  0)  is  the  complex  function  defined  in  Equation  [34].  The 
e  (0),  and  are 


functions  f{k,a,z). 


[67] 


V  0 

•  f(k,a,z)  =  kh  cosh  ka  cosh  kz 

e{0)  =  cos^  0 

g(k)  =  cosh  kh  sinh  kh 

The  functions  w)  and  /gCi?,  w)  are  defined  in  Appendix  B.  The  functions  N ^  and  Aj 

needed  for  w)  are 


/V,  =  AA  -* /I4Z4  -  A3Z4  - /1^Z3] 

Aj  =  Kh  U3Z3  -  /I4Z4  +  /I3Z4  +  /I4Z3] 

’  d 

In  the  supercritical  range  of  Froude  number  it  is  necessary  to  determine  the  Cauchy  principal 
value  at  the  singularity  which  occurs  when  F  cos  (9  =  1.  Explicit  forms  of  K (f,)  are  the 
same  as  those  for  h'y^irf).  Equations  [40]  through  [44],  if  the  functions  g{e)  and  /(/fc^)  are 
•  given  the  values 

•  • 

(AqA)^  cos^  0  cosh  kf.a  cosh 

(")  =  — - — ^ - V-  [69] 

sinh  AqA  cosh  k^h  [F^  cos^  0  -  sech^  /fc^A] 


/(^o) 


{k^h)'^  cos^  0  cosh  kf^a  cosh  k^z 
sinh  AqA  cosh  k^h  sin  0 


[70] 


The  derivatives  of  are  obtained  from  the  derivatives  of  the  kernel* functions.  The 
three  derivatives  of  (,,)  in  Equation  [65]  are 


dK, 


dx 


^  L  J 


[71] 


dK, 


V 


I  3 


dKj^ 

l^^enh-a  z.2nh^al 


.[72] 


[73] 


where  R^  and  R^  are  as  defined  in  Equation  [47].  The  derivatives  of  the  other  kernel 
functions  are  obtained  from  the  derivatives  of  (r,)  in  Equation  [66].  The  a:-derivative  is 
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the  real  part  of 


dK, 


=  -•/{'/.  «>+)  -J  (t),  w_) 


[74] 


where  e  (6)  is  the  product  of  cos  d  and  its  former  value.  The  y-derivative  is 


-  J  (r),  u>_)  [75] 


where  e  (d)  is  the  produot  of  sin  0  and  its  former  value.  In  both  cases  the  new  value  of 
f(k,a,s)  is,  ^ 

m 

f{k,a,z)  =  {kh^^  cos\\Jca  cosh  kz  [76] 

Expressions  for  the  functions  ./j  (»?,  w)  and  ^3(77,  w)  may  be  found  in  Appendix  B.  The  func¬ 
tions  N  j  and  needed  for  J2(Vi  “r® 


Nl=-2(hhf  [A^Z^  +  A^Z^] 
2(A'A)M/l3Z3-/l,Z^] 

The  3-derivative  of  Up  (tf)  is 

dK  D 

^  '  /('?.  +  l(v,  wj 

The  function  e  (6)  is  unchanged  but  f(k,a,z)  becomes 

f(k,a,z)  =  (kh)^  cosh  ka  sinh  kz 
The  functions  /Vj  and  needed  for  u)  are 

A/j  =-2(0)2  [^3^2*+  A^Zj] 

A/2=  2(0)2  U3Z, 


•  [77] 


[78] 


[79] 


[80] 
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STRENGTHS^  OF  TH&  SINGULARITIES  USED  TO  REPRESENT  A  HYDROFOIL 

The  strength  of  vorUcity  along  the  span  of  a  hydrofoil  is  a  function  of  the  chord  and  ^ 
lift  developed  at  any  section  and  the  free-stream.velocity.  If  the  vorticity  is  assumed  to  have 
an  elliptic  distribution,  the  vorticity  per  unit  length  V  (r,)  in  the  equations  for  has  the  form 


2L  !  ^  o  cC ,  U  j  2 

r(,,)=_^  l/i -1- =  Z  __L.  i/i  _:l 

y  n  h  y  ^  .2 


[81] 


In  this  expression  L  is  the  total  lift  and  is  the  average  lift  coefficient  of  the  foil. 

The  strength  per  unit  length  of  the  horizontal  line  doublet,  used  to  represent  the  dis- 

p  acement  effect  of  the  hydrofoil,  may  be  written  in  terms  of  a  shape  parameter  A  as  follows 

• 

[82] 

oA. 


where  t{r,)  is  the  maximum  width  of  the  hydrofoil  section.  If  the  hydrofoil  had  a  circular 
cross-section,  A  would  be  unity  and  would  be  the  usual  strength  of  a  line  doublet 

Lock  has  compyted  A  for  a  number  of  different  geometric  cross  sections.  For  typical  hydro¬ 
foil  sections  A  is  approximately 


9  A  =  0.42  (1  +  c/t)  [83] 

•  *9  m 

where  c/t  is  the  length-width  ratio  of  the  section. 
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appendix  a 

INTEGRAL  RELATIONS  USED  IN  JHE  POTENTIAL  FiTncTIONS 

The  integral  relations  used  for  the  potential  functions  in  the  foregoing  analysis  are 
discussed  in  'Vu’s  report^  and  arc  summarized  here.  If  the  relations  are  known  for  one  of  the 
singularities,  the  others  may  be  obtained  by  differentiation  or  integration.  If  the  basic  inte¬ 
gral  relation  is  given  for  the  source,  the  doublet  relations  are  obtained  by  differentiating  this 
function  with  respect  to  one  of  the  coordinates.  The  integral  relation  for  the  semi-infinite 
vortex  sheet  is  obtained  by  integrating  the  expression  for  the  vertical  doublet  from  the  posi¬ 
tion  of  the  hydrofoil  downstream  to  infinity. 

As  the  potential  of  a  source  is  proportional  to  the  reciprocal  of  the  dfstance,  the  first 
integral  relation  is  • 


/ar^  + 


where  the  zeroth-order  Bessel  function  {k  ^  y^)  may  be  defined  by  the  integral 

*  2  ''/2 

<^0  +  y^)  =  ~  J  cos  (kx  cos  d)  cos  (ky  sin  9)  dO 

0 

1 

~  ~  J"  fcos  kWj^  +  cos  k‘W_]  dO  [2A] 


where 


w-V  =  X  cos  6?  +  y  sin  0 


=  X  cos  0  —  y  sin  0 


Equation  [2A]  may  also  be  written  as  the  real  part  of 


Jq  {k  x/a:^ 


-  1  ^  •  I  2 

y^)  -Ke'-  J  («■*“■-  * 


In  the  limit  as  approaches  infinity,  this  integral  becomes 


2  f  sin  u  * 

—  cos  ky  I  -  du  =  cos  A?/ 

n  J  u 


Therefore,  the  two  integral  relations  for  the  semi-infinite  line  distribution  of  doublets  are  the 
real  parts  of 


OO 

-  =Re  sgn  z  f  dk 

+  2^  -I 


n/ 2 

sgn  z  r  do 


+  2^  /  2  2  2 


/*  I-  I  I  ikw  , 

- -  /  2  e  +  dk  [12A] 

COS  0  J 


When  Boundary  Condition  II  is  applied  in  the  case  of  the  semi-infinite  line  distribution 
of  doublets,  the  kernel  function  (^)  from  Equation  [19]  is 


,  n/2 

h  r  do 


r  /  \  ti  ^  f  C  ikw  ,, 

-nf  2  0 


For  subcritical  Froude  number 


F{k,a,z)  =  agn  {z-a)  _  ^-k  u  ^  a)  ^  ^-kh  sinh  ka  cosh  kz 

sinh  kh 

If  the  hyperbolic  functions  are  written  in  their  exponential  forms 


,a,2)  =  ^  sgn  (2nA  +  2-0)  "A  +  z  -  «l  _  ^~’e~ 

M  ~  n 


k  {2  nh  —  z  +  a) 


«>  «> 

-  ^  'g~*  (2  nA  +  z  +  o)  ^  \  'g-k  (2  n 

«  =  0 


h  —  z—  a) 
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If  the  sign  of  n  is  changed  in  the  second  and  fourth  sums,  the  terms  combine  to 


^  OO 

F{k,a,z)=  ^  sgn  (2nA  +  2-a)  ^  sgn  (2nA+  2  +  a)  c"*! 2  "A  +  z  +  o| 


If  the  integral  relation  of  Equation  [12A]  is  used 


n  ^  —  OO 


2nh  z  —  a 


(r,-yf.{2nh^z-a)^  ^^2  ^  /(2„A+ 2-a)2 


2  nA  +  2  +  a 


('/~y)^  +  (2nA  +  2  +  a)2  /  2  ^  /  Ti  7  73 

'  V*  +  (T/-y)  +  (2nA +2  +  a)2 


In  the  supercritical  range  of  h’roude  number,  the  last  term  in  the  expression  for  F  (k,a,z)  is 
-cosh  kh  instead  of  sinh  kh;  see  Equations  [25]  and  [26].  Then  in  the  kernel  function 
('/)  terms  in  the  infinite  sums  alternate  in  sign. 


APPENDIX  B 

EVALUATION  OF  THE  INTEGRALS  BY  CONTOUR  INTEGRATION 


The  wavemaking  potentials  derived  in  this  report  are  obtained  from  the  real  or  imagi 
nary  part  of  an  integral  of  the  type 

,  =  </(n.  +  ■/(t?,  MiJ  -  t  [/(»?,  +  /(,,,  «,_)]  [IB] 


Whereas  the  kernel  functions  of  the  potentials  K  (tf)  and  their  2-derivatives  are  of  the  type 
/  (r,,  w),  the  and  y-derivati  ves  are  of  the  type  J  (rj,  tc).  In  this  expression  J  (r,,  u)  and 
/(rf,  w)  are  the  real  and  imaginary  parts  of  the  complex  integral 

n/2 

J (t),  w)  -  i  I  (rf,  w)  =  J  e{0)  M  {w,  6)  de  ®  [2B] 

0 

where  0)  is  tHe  principal  value  of  the  integral 


”  J  ^ ,  V  r_  tan ! 


0  ff  (^)  \F^  cos^  0 


tanh  kh' 
kh 


[3B] 


In  these  equations  and  w_  are  as  defined  in  Equation  [32],  and  w  is  written  for  either 
or  u_.  The  functions  e{0i)  and  depend  upon  the  type  of  singularity,  and  g  {k)  is 


g(k)  =,sinh  kh  cosh  kh  [4B] 

The  integrals  are  not  defined  for  F  -  1. 

The  integrals  in  Equations  [2B]  and  [3B]  have  singularities  in  the  integrands  and 
must  be  evaluated  by  means  of  a  contour  integration.  In  the  integral  defining  M{w,  0)  there 
are  no  singularities  in  the  quotient  f{k,a,z)  'g{k).  There  is  a  simple  pole  in  the  integrand, 
however,  when  k  =  Aq,  where  is  defined  by  the  equation 

•  • 

m 

kgh  F^  cos^  0  =  tanh  k^h  [5B] 

In  the  complex  A:-plane  there  are  infinitely  many  singularities  on  the  imaginary  axis,  and  the 
contour  of  integration  must  be  chosen  in  such  a  way  that  this  axis  is  excluded  and  the  inte¬ 
grand  remains  finite.  If  w  is  positive.  Contour  I  of  Figure  2,  which  is  in  the  upper  part  of 


Contour  I 


w>'0 


Contour  II 


M<0 


Figure  2  Contours  in  the  Complex  A-Plane  for  Kvaluating  the  Integrals 

A 

the  plane,  may  be  used.  If  w  is  negative,  Contour  11  in  the  lower  part  of  the  plane  may  be 
used.  The  signs  of  and  fo_  for  different  ranges  of  x  and  (9  are  shown  in  Table  1  for  posi¬ 
tive  values  of  T  -  y.  For  negative  values  of  r?  -  y,  and  w_  are  interchanged.  In  this  table 


tan0,= 


[6B] 


TABLE  1 

Signs  of  tc^  and  u.'_  for  Different  Ranges  of  x  and  0 


^  0 

0 

<  0  < 

n/2 

Contour  I 

X  >  0  •< 

«J_  >  0 

0 

<  0  < 

77/2  - 

Contour  I 

L'"_  <  0 

77/2  -  flj 

<  0  < 

77/2 

Contour  II 

• 

rw_  <  0 

0 

<  0  < 

77/2 

Contour  II 

X  <  0  ^ 

1  «; .  <  0 

1 

0 

<0  < 

77/2  - 

Contour  II 

>  0 

n/2  -  (9j 

<0  < 

77/2 

» 

• 

Contour  I 

Figure  2  shows  the  pole^within  Contour  I  at  Atq  +  ia,  where  a  approaches  the  limit  zero. 
Therefore  in  the  integration  about  Contour  I,  there  is  a  residue  term  which  provides  wave- 
making;  but  there  is  no  residue  or  wavemaking  associated  with  Contour  II.  For  large  positive 
values  of  x,  is  very  small  and  and  u_  are  positive  over  most  of  the  range  of  0.  For 
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smallTiegative  values  of  ar,  there  is  wavemaking  ahead  of  the  hydrofoil  which  decreases  as  x 
increases  and  vanishes  as  %  approaches  minus  infinity.  Therefore,  Boundary  Condition  III  is 
satisfied.  If  the  pole  were  exactly  on  the  axis  and  shared  equally  by  both  contours,  there 
would  be  as  much  wavemaking  ahead  of  the  hydrofoil  as  downstream  from  it.  This  clearly 
contradicts  Boundary  Condition  III,  and  Figure  2  gives  the  correct  location  of  the  pole.  This 
analysis  is  equivalent  to  the  use  of  a  fictitious  viscosity  to  locate  the  singularity. 

When  w  is  positive,  the  integration  around  Contour  I  gives 


a  a  a 

j^[  ]dk  =  -J  [  ]ke‘^dO^  J  [  ]  d  \k\  +  2ni  Resik^)  [73] 


r  9 

In  the  limit,  when  A  and  B  become  infinite,  the  integral  over  the  arc  AB  vanishes.  Along  OB 

* 

Ijfcl  +  [8B] 

When  w  is  negative,  the  integration  around  Contour  II  gives 


a  n 

I  ]  dk  =  J  [  ]  (I- i)  dK 


If  the  term  for  the  integration  along  OB  is  given  the  subscript  2  and  the  term  for  the  residue, 
subscript  3,  then  M  is  given  by  the  sum  of  M2  and  M where 


hr  f  (K  {li-isgnw),  a,z)  e 

M2(uj,0)  =  —  I  - 

0  y  (A’ (1  +  t  sgn  to)  )  cos^  0 


sgnto),  a,z)  e  e ‘Aw (1  +  isgn  w)  (iA' 


\f^  cos 2  p  tanh  (FA  (l  +  tsgnw))1 
L  h'h  (1  + 1  sen  w)  J 


-  i  sgn  tc 


00 

2A  r 

n  J  q  (K  (\  +  ii 


/(A' (1  +  tsgn  M>),  a,z)  e  dK 


J  y  (A"  (1 +  tsgn  tc)) 

0  - 777 -  lA'A  (1  +  tsgntc)  cos2  0  -  tanh  (Kh  (1  +tsgntc))] 


W3  (tc,  0)  = 


t  (1  +  sgn  tc)  f  (Aq,  a,  z)  e 
g  (A(j)  IF^  cos 2  0  -  sech2 
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To  obtain  the  real  and  imaginary  parts  of  6),  let 


A/,  +  t  sgn  w  N 2  =  f  {K  (l>isgnM)),  a, 2) 


.  g  {K  (l  +  isgnwj))  -i  o 

t/j  +  t  sgn  w  £>2  =  - — -  [Kh  (1  +  isgntc)  cos*'  6  -  tanh  {Kh  (1  +  isgntr))] 


[13B] 


The  real  and  imaginary  parts  of  ^2(1^)  are 


2/1 


J2iri,  w)  ^  J  e(d)dO  J  IC,  sgn  «)  sin  /(■«  +  Cj  cos  A'u;|  rfA"  [14B] 


TT  f  A  ao 

•  2h  r  r 

’")  =  ■—  J  e(0)d0  J  \C  ^  sgn  w  cos  Kw  -  C2  sin  Kw\  dK  [15B] 


where 

•  A/,«,  +  /V2f>2  ^2  ^1-^1  ^2 

C^ - - - - - ;  C2  =  -  [16B] 

The  values  of  and  niust  be  determined  for  the  particular  singularity  under  consideration. 
These  are  functions  of  the  following  quantities: 


/Ij  =  sinh  Ka  cos  Ka  .  ^2  = 


^43  =  cosh  Ka  cos  Ka 


Z  J  =  sinh  Kz  cos  Kz 
Z3  =  cosh  Kz  cos  Kz 


=  sinh  Ka  sin  Ka 


7^2  =  cosh  Kz  sin  Kz  1 
=  sinh  Kz  sin  Kz  ) 


Then  the  expressions  for  //j  and  A/j 


^  ^  ("1  -"2)  -  "3  +  1] 


=  2A'A  0(1!  2)- 
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where 


Wj  -  sinh  2A’4  cos  2Kh 
W3  •=  cosh  2Kh  cos  2Kh 


^  //j  “  cosh  2Kh  sin  2Kh 
=•  sinh  2KJi  sin  2Kh 


As  K  becomes  small  and  in  Equation  [16B]  are  of  magnitude  {Kh)’',  where  n  is  equal 
to  or  greater  than  unity  for  ail  the  singularities  considered.  In  the  denominator  Z?  j  and  Z? 
have  magnitudes  {Kh)^  and  Kh  ^  cos^  6  -  1],  respectively.  Therefore,  Cj  and  are  de¬ 
fined  for  alWalues  of  K  and  6  in  the,subcritical  range  of  Froude  number.  In  the  supercritical 
range  of  Froude  number,  F  cos  0  may  also  be  greater  than  unity  and  becomes  equal  to  unity  « 
when  0  =  0^,  defined  by- 

cos  Oq  = 1/F  [21B] 

When  the  value  of  the  integrand  is  investigated  as  6  approaches  6^  from  both  sides  of  the 

singularity,  the  Cauchy  principal  value  is  found  to  exist  and  the  integrals  in  Equations  [14B] 
and  [15B]  are  defined. 

The  real  and  imaginary  parts  of  the  functions  defining  for  the  subcritical  range 

of  Froude  number  are  * 


/T  /  2 

r  «  (<?) /(Z:o,a,z)  (1+ sgnw)  sin  kr.w  k^de 

■^3  = -A  /  — - - ; - - - ° - 

J  sinh  k^h  cosh  k^h  cos^  e  -  sech^  /t^A] 


/  2 

,  ,  ,  r  ^  (®)  /  (1  +  Sgn  w)  cos  AoW  A-de 

— - - - ^ - 2 - 2 -  (23B, 

0  sin^  k^h  cosh  k^h  [F^  cos/  0  -  sech^  k^h] 

As  long  as  the  Froude  number  is  less  than  unity,  the  smallest  value  k^  can  assume  is  always 
greater  than  zero  and  these  expressions  are  defined  for  all  values  of  k^  and  0.“ 

In  the  supercritical  range  of  Froude  number,  the  expressions  for  J,{t,,w)  and  Z,(n,w) 
are  the  same  except  for  the  lower  limit  of  the  integral  which  is  0^  instead  of  zqro.  As  k^  is 
zero  when  0  =  0^,  these  expressions  appear  to  have  singularities  for  these  values  of  the  pa¬ 
rameters.  The  difficulty  is  avoided,  however,  if  the  integration  is  evaluated  over  k.  instead 
of  0.  If  Equation  [5B]  is  differentiated 


F2cos2  0_sech2AoA  2F2cos0sin0 
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then 


j'(ri,w)  =  -  -AL  f  ‘  (1  +  sgnw)  sin^oM  k^dk^ 

2F^  "I  sinh  AtqA  cosh  AqA  cos  0  sin  0 


[25B] 


^3  ('?»«’) 


J  *  /(Ao.«.2)  (1  +  sgnM)  cos  AqM 


sinh  AqA  cosh  ijjA  cos  0  sin  0 


[26B] 


The  integrands  of  these  integrals  have  no  singularities  and  converge  to  zero  as  k 
infinity  and  0  approaches  n/2.  ' 

The  derivatives  of  C  (,)  with  respect  to  y,  and  2  are 


approaches 


A 


dL 

dx 


=  f(ri,  W^)+  t 


[27B]. 


where  e(9,  ,,  ^ 


A 


dL 

—  =-  /(,,  +  /(,,  „,j  +  i  [j  ] 

^  •  • 


[28B] 


where  e{6)  is  replaced  by  the  product  of  sin  0  and  its  former  value, 
becomes  the  product' of  k  and  its  former  value. 


In  both  cases  f{k,a,z) 


dL 

dB  ^  <>'('?.  «+)  +  J{v,  wj)  -  i  [/(^,  +  i^r,,  «;_)  ]  [29B] 


where  e(0)  is  unchanged  but  f{k,a,B)  is 


replaced  by  the 


2-derivative  of  its  former  value. 
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